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We study a magnetic impurity embedded in a correlated electron system using the density-matrix 
renormalization group method. The correlated electron system is described by the one-dimensional 
Hubbard model. At half filling, we confirm that the binding energy of the singlet bound state 
increases exponentially in the weak-coupling regime and decreases inversely proportional to the cor- 
relation in the strong-coupling regime. The spin-spin correlation shows an exponential decay with 
distance from the impurity site. The correlation length becomes smaller with increasing the corre- 
lation strength. We find discontinuous reduction of the binding energy and of spin-spin correlations 
with hole doping. The binding energy is reduced by hole doping; however, it remains of the same 
order of magnitude as for the half-filled case. 
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I. INTRODUCTION 

Although more than forty years have passed since the 
discovery of the Kondo effect, it is still one of the most 
interesting topics in condensed matter physics; it lies at 
the heart of understanding strongly correlated electron 
systems^ The Kondo effect, which leads to the quenching 
of an impurity spin, forms the basis of the physics of a 
single magnetic impurity embedded in a metal. In order 
to understand the Kondo effect, the Anderson model 2 
has been applied with great success. In the theoretical 
studies, one generally assumes an impurity level to be 
embedded in a noninteracting conduction band. 

In the past a system of magnetic ions coupled to 
(strongly) correlated conduction electrons has attracted 
considerable interest in connection with heavy-fermion 
behavior, namely Nd2- a; Ce a ;Cu044 This raised the ques- 
tion whether correlations among conduction electrons af- 
fect substantially the expected formation of heavy quasi- 
particles^ So far a number of authors have studied mod- 
els of a single magnetic impurity embedded in a host 
of correlated conduction electrons. Thereby perturba- 
tion theory and other approximation schemes were ap- 
plied ,£i£iLSi2 For example, it was shown that the Kondo 
scale can increase exponentially in the weak-coupling 
regime with increasing interaction of the conduction elec- 
trons!^ But a quantitative theory is still missing. More- 
over the case of strongly correlated conduction electrons 
with band filling slightly less than one-half (hole doping) 
is still an open problem. 

In this paper we study a single magnetic impurity 
coupled to a correlated electron system. The latter 
is assumed to be one dimensional (ID) and described 
by a Hubbard Hamiltonian. Using the density-matrix 
renormalization group (DMRG) method, we calculate 
the binding energy of the impurity-induced bound state 
and spin-spin correlation functions between the impurity 
and the correlated electrons in the thermodynamic limit. 
Special attention is paid to the case of a nearly half-filled 
conduction band with repulsive electron-electron interac- 
tions. For a ID correlated host, there have been numer- 



ical study for similar models^ as well as an analytical 
study for a integrable models We hope that the present 
investigation will contribute to better insights. 

This paper is organized as follows. In Sec. II, we in- 
troduce our model, i.e., a magnetic impurity coupled to 
a Hubbard chain. In Sec. Ill, we give some numerical de- 
tails of the DMRG method applied here. In Sec. IV, we 
first present calculated results for the binding energy and 
spin-spin correlation functions at half filling, and discuss 
the effect of the host-band correlations on the Kondo 
physics. Then, we consider the evolution of the same 
quantities with hole doping. Sec. IV contains a summary 
of the results and the discussions. 



II. MODEL 

We study a magnetic impurity coupled to a ID corre- 
lated electron system. The Hamiltonian consists of three 
terms 

H = H c + Hf + H c f. (1) 

The first term H c represents ID correlated electrons. 
Here we describe them by a Hubbard Hamiltonian, 

He = *l>j+ 1(7 c iCT + h.c.) + Uj2 n H n ih ( 2 ) 

i,a i 

where c\ a (ci a ) is the creation (annihilation) operator of 
an electron with spin a (=T, I) at site i, and m a = c\ a Ci,j 
is the number operator. Furthermore t is the hopping 
integral between neighboring sites and U is the onsite 
Coulomb interaction. The second term Hf is the orbital 
energy of the magnetic impurity site. We assume that the 
impurity contains one orbital, e.g., 4/ and the Coulomb 
repulsion on the orbital Uf is infinite. Since double oc- 
cupancies are excluded, i.e., the / orbital is either empty 
or singly occupied, the impurity site is given by 

H f = s f J2PM (3) 
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FIG. 1: Lattice structure of the system. Open and solid 
circles represent correlated electron system and impurity site, 
respectively. The bottom numbers i denote site index of cor- 
related electron system and |i| corresponds to a distance be- 
tween site i and the impurity site. 

with fl — — flfv) and Ef < 0. For convenience, 
we define r = —Ef/U(> 0) and label an electron on the 
impurity site as "/ electron" . The third term H c f in- 
volves the interaction between the impurity site and the 
correlated electron system. The interaction is assumed 
to be local and described by a hybridization like in the 
Anderson model, i.e., the impurity site is hybridized with 
a single site (denoted as site 0) of the correlated electron 
system. Thus, 

H cf = vY,{clJ« + fl^l (4) 

a 

where V is the hopping integral between the impurity 
site and site 0. The lattice structure is shown in Fig. [1] 
We will work in units where t = 1 and take as typical 
values |e/| =2 — 3 and V = 0.1 — 0.2, throughout. 

III. METHOD 

We employ the DMRG method, which is one of the 
most powerful numerical techniques for studying quan- 
tum lattice many-body systems including quantum im- 
purity systems^ With the DMRG method, we can ob- 
tain ground-state and low-lying excited-states energies, 
as well as expectation values of physical quantities quite 
accurately for very large finite-size systems. 

In order to carry out our calculations, we consider N 
(— Ni + Ni) electrons (N: even) in a system consisting 
of a chain of Lsites correlated electron system {L: odd) 
and a single impurity site. The electron density is defined 
as n = N/ (L + 1). Note that the number of lattice sites 
must be taken as L + 1 = 41 — 2 with I (> 1) being an 
integer to maintain the total spin of the ground state as 
5 = 0. If one chooses it as L + 1 = 4Z, the singlet and 
triplet states are degenerate. We now apply the open- 
end boundary conditions to the ID correlated electron 
system and assume that the impurity site is hybridized 
with the central site of the ID open chain. The latter 
corresponds to site 0, and sites i and —i are equivalent. 
In this paper we restrict ourselves to the half-filled and 
hole-doped cases (N < L + 1). 

Regarding quantum impurity problems, it is generally 
complicated for finite-size calculations to obtain accurate 
results in the thermodynamic limit L — > oo because of 



finite-size effects. In our calculations, the most prob- 
lematic finite-size effects arc Fricdcl oscillations due to 
the open ends of the Hubbard chain. Mostly, the energy 
scale of the Kondo physics is exponentially small; never- 
theless, Friedel oscillations can persist even at the center 
of the chain as they decay as a power-law from the edge 
sites. Therefore, we study several long chains with sites 
L + 1 = 62, 126, 190, 254, 318, 382, 446, and 510, and 
then perform the finite-size-scaling analysis based on the 
size-dependent quantities. All DMRG results in this pa- 
per arc extrapolated to the thermodynamic limit L — ► oo. 
For precise calculations we keep up to m ss 5000 density- 
matrix eigenstates in the DMRG procedure. In this way, 
the maximum error in the ground-state energy is below 
10" 8 - 10" 7 . 



IV. RESULTS 
A. System at half filling (n = 1) 

1. Binding energy 

We first study the binding energy between the / elec- 
tron and the correlated electrons. It corresponds to an 
energy gain due to the formation of a Kondo (or local) 
singlet bound state. Hence, the binding energy is given 
by an energy difference between the first triplet excited 
state and the singlet ground state, 

A B = lim A B (£), (5) 

L ►oc 

with 

A B (L) = E (L, iV T + 1, JV; - 1) - E (L, N h N t ) (6) 

where Eq(L, Nf,N±) is the ground state energy in a sys- 
tem of L + 1 sites with N-\ up-spin and N± down-spin 
electrons. Note that, at half filling, the system is insulat- 
ing for finite U . The bound state therefore may be from 
as a local singlet rather than the Kondo singlet. Here and 
in the following we will speak of a Kondo singlet only if 
it involves more than the central site of the correlated 
electrons. 

In Fig.[5]Ja), we show the DMRG results of the binding 
energy Ab as a function of the Coulomb interaction U 
for various parameter sets. In total, the results for the 
different parameter sets are qualitatively the same; as U 
increases, the binding energy rises rapidly for small U, 
reaches a maximum around U « 4 and decreases grad- 
ually for large U. This behavior is similar to the de- 
pendence of the effective Heisenberg interaction on the 
Coulomb interaction in the half-filled Hubbard models 
Accordingly, the DMRG results show that for large val- 
ues of U the binding energy is approximately propor- 
tional to the effective exchange coupling J c f, between 
the impurity and site 0^ If we assume that the effec- 
tive exchange coupling results from second-order pertur- 
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FIG. 2: (a) Binding energy A B for e } = —3, V = 0.2 
(circles), sj = —3, V = 0.1 (triangles), and £/ = —2, 
V = 0.1 (squares), (b) Magnified view of small U region 
for Ef = —3, V = 0.2. The data is fitted by a function 
A B = \ZaUexp{-/3/U) with a ~ 2.5 x 10 -4 and /3 ~ 0.4. 
(c) Ab for e/ = —3, V — 0.2 in the strong-coupling regime 
([/ > i). The data is fitted by a function Ab — u 2 £ with 

7~3.1xl0~ 4 . 



bation, i.e., J c / = ^rg~) , we can explain why the re- 
sults for V = 0.2 are about four times larger than those 
for V = 0.1. This estimation of the effective exchange 
coupling is also consistent with a slight decrease of the 
binding energy with increasing \sf\. 

Let us now consider the behavior in the limiting cases 
for weak and strong interaction strengths. A mag- 
nified view of the weak-coupling regime (U < t) for 
e f = -3,V = 0.2 is given in Fig. E^b). When U = 0, 
the system is metallic and essentially equivalent to the 
single-impurity Anderson model (SIAM) in the Kondo 
limit (Uf/V = oo ) but asymmetric case (e ^ — C///2). 
The orbital energy of the impurity site is lower than the 
Fermi energy of the conduction band, so that the oc- 
cupation number of the impurity site is always 1. The 
exchange interaction J c f is estimated to be the order of 
V 2 /sf and, therefore, the binding energy is expected to 
be very small but finite. We estimate it to be roughly 
Ab — 10~ 7 — 10~ 6 . This value is compatible with the 
Kondo temperature Tk in the asymmetric SIAMJ^ The 
introduction of a finite Coulomb interaction makes the 
system insulating. With increasing U, Ab increases grad- 
ually when U/t < 0.2 and rapidly for U/t > 0.2. There 
is a crossover from the Kondo singlet to a local singlet 
around U/t — 0.2. Assuming an exponential behavior of 
Ab with U leads to a good fitting of the DMRG data, i.e., 
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FIG. 3: Spin-spin correlation functions (a) (S/ • so) and (b) 
{Sf ■ si) as a function of the Coulomb interaction U for £f = 
—3, V — 0.2 (circles), Sf = —3, V — 0.1 (triangles), and 
£/ = —2, V — 0.1 (squares). Inset: semilogarithmic plots of 
the magnitude of the spin-spin correlation functions. 

A B = VaUexp(-/3/U) with a ~ 2.5 x 10~ 4 and (3 ~ 0.4. 
Furthermore, Ab increases almost linearly in the regime 
U/t = 0.2 — 2. We thus find that the binding energy of 
the local singlet can be a few orders of magnitude larger 
than that of the Kondo singlet. 

The DMRG results for the strong-coupling regime 
(U ^> t) with Ef — —3 and V — 0.2 are plotted in 
Fig. Etc). In this regime, the electrons are strongly lo- 
calized at each site. Therefore, the system ([1} can be 
reduced to the Heisenberg model with Hamiltonian 

H cS — J Si ■ s i+ i + JcfSf ■ s (7) 

i 

with J = The DMRG data can be fitted quite well 
by a function Ab = with 7 ~ 3.1 x 10~ 4 . Despite 

the strong localization of the electrons the binding energy 
is two orders of magnitude smaller than the cf exchange 
coupling. This is because for n = 1 a spin-density-wave 
(SDW) is forming in the chain for any value of U (> 0) 
which makes the formation of the local singlet state more 
difficult. We note that the behavior of the binding energy 
for finite U is essentially the same as that of the Neel 
temperature in the half-filled Hubbard model^ 

2. Spin-spin correlations 

In the Kondo problem, the spin degrees of freedom 
around the impurity play an essential role. Therefore, 
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we investigate spin-spin correlations between the / elec- 
tron and the correlated electrons. The correlated system 
is now described by the lattice model [Eq.(j2}], so that 
we are allowed to study the distance r dependence of 
correlation functions, like (Sf • s r ). 

Let us first derive the spin-spin correlations between 
the spin on the impurity site and on the central site i = 0, 
i.e., (Sf ■ So). The DMRG results for various parameter 
sets are shown in Figj3^a) as function of the Coulomb 
interaction U. Since the antiferromagnetic correlation 
is derived from the cf exchange interaction, (S f ■ Sq) is 
negative for all parameter sets and Coulomb interaction 
strengths. The absolute value of (Sf ■ Sq) increases with 
increasing V and with decreasing |£p|, as expected from 
the behavior of the binding energy. However, the influ- 
ence of |ep| is rather smaller. In the limit U — * 0, (Sf ■ So) 
is antiferromagnetic but the magnitude is very small due 
to strong charge fluctuations, when the system is metallic 
[see inset of FiglH^a)]. It reflects the small binding energy 
around U — 0. The magnitude of (Sf ■ sq) increases with 
increasing U and reaches its maximum value as U — > oo, 
which means that one electron is localized on each site in 
that limit. 

We consider next spin-spin correlations between a spin 
on the impurity site and on the next-nearest-neighbor 
site i = 1, i.e., (S f • s x ). In FigG^b), the DMRG results 
for (Sf • Si) are shown as a function of the Coulomb in- 
teraction U for various parameter sets. One expects fer- 
romagnetic correlations from the effective Hamiltonian 
(J?]) for finite values of U, and indeed (Sf • si) has posi- 
tive sign for all the parameter sets and Coulomb interac- 
tion strengths. Note that the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction induces ferromagnetic corre- 
lations, as substitute for the spin-spin interaction (J7J, in 
the weak-coupling (U ~ 0) and metallic regimes. How- 
ever, it is difficult to separate the contribution from 
RKKY and the interaction ((7J . The Coulomb interaction 
dependence of (Sf ■ si) is similar to that of (Sf ■ So). For 
the same parameter sets, the value of (Sf ■ si) is found 
to be slightly smaller than that of |(S/ • so)|. This indi- 
cates a slow decay of the spin-spin correlation (S/ • s r ) 
with distance r. It implies that the spin of the / elec- 
tron is hardly screened by the spin on site 0. In addition, 
the influence of V on the spin-spin correlations is rather 
small. Note that the binding energy depends strongly on 
the hybridization V . 

Let us now consider the distance dependence of the 
spin-spin correlation functions. In Fig. 01a), we plot 
the DMRG results for (S function of distance 

r (= We choose three Coulomb interactions; (i) 

U = 0.5 in the Kondo-singlet regime, (ii) U — 200 in 
the limit of the local singlet regime, and (iii) U = 4 in 
the intermediate regime where a maximal binding en- 
ergy is obtained. The results for different distances are 
extrapolated to the thermodynamic limit L — ► oo . We 
find that (Sf ■ s r ) decays slowly and the sign changes al- 
ternately with r, i.e., (Sf ■ s r ) has a positive (negative) 
sign for odd (even) r, denoted by solid (empty) sym- 
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FIG. 4: (a) Spin-spin correlation functions (Sf ■ s r ) as a func- 
tion of the distance r for U = 0.5 (triangles), 4 (squares), and 
200 (circles), (b) Semilogarithmic plot of the magnitude of 
the spin-spin correlation functions. The data is fitted by a 
function (Sf ■ s r ) ~ exp(-|) with £ = 3184, 508, 400 for 
U = 0.5, 10, 200, respectively. 

bols in Fig. |4ja). The interaction flTJ) and/or the RKKY 
interactions cause ferromagnetic (antiferromagnetic) cor- 
relations between the spin of the / electron and that of 
the odd (even) site r. The absolute value of (Sf ■ s$) in- 
creases with increasing U because larger Coulomb inter- 
actions stabilize the 2fcF-SDW oscillation which accom- 
panies charge localization. 

Since the system is in a spin-gapped ground state, 
an exponential decay of the spin-spin correlation with 
distance must be expected. In Fig. |4jb) we present a 
semilogarithmic plot of (Sf ■ s r ) as a function of distance 
r. For a convenient comparison, we have normalized 
(S f ■ s r ) with respect to its value at r = 0. The results 
can be fitted with a function exp(— |) and thus the ex- 
ponential decay of the correlation functions is confirmed 
for all values of U. The correlation lengths are estimated 
as f = 3184, 508, 400 for U = 0.5, 4, 200, respectively. 
They seem to be much longer than those of other stan- 
dard spin-gapped systems, e.g., £ = 3.19 in the two-leg 
isotropic Heisenberg system. The large values of £ reflect 
exponentially small binding energies. They also mean 
that spin-polarized electrons are widely spread around 
the impurity site, i.e., the Kondo screening effect is quite 
weak. Furthermore, we note that the correlation func- 
tions decay rapidly around r ~ 0. The decay rate is 
dependent of the magnitude of the cf exchange interac- 
tion. 
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B. Less than half filling (n < 1) 

We are also interested in doped systems, which are 
metallic even if U > 0. We thus investigate the prop- 
erties of the model ([1]) with £/ = —3 for various hole 
concentrations n = 1 — N^/L, where is the number 
of doped holes (N^ > 0). For this choice of the occu- 
pation number of the impurity site is near unity because 
the Fermi level lies well above In the strong cou- 
pling limit (U ^> i), doubly occupied sites are excluded 
and therefore we can derive an effective model ([TJ by 
applying degenerate perturbation theory— The effective 
Hamiltonian is written as 



[X10" 4 ] 



H = H t +Hj + H p + H K + H'. 



(8) 



Here H t is the kinetic-energy term of the conduction elec- 
trons 

t ( V 2 {2 + 2r + r 2 ) 



2e 2 {l + r) 2 



S i0 , (9) 



with c| CT = c\ a (l — riiu). Furthermore Hj is a spin- 
coupling term between the conduction electrons, which 
is of the Heisenberg type 



Hj — Ji Si • Sj+i, 



2t 2 , 

Ji = 1 



V 



U V £f(U-s f ) 



>i.Q 



(10) 



The sum of these two terms defines the ID correlated 
electron system. It is essentially equivalent to a t—J 
model except for small modifications around site due 
to the impurity. The term H p corresponds to one-particle 
potential around the impurity site, which is given by 



7]V 2 



Hp = 25,(1 + r) (1 " U0) + 7m + r) 2 X> " 



-/ 



i=±l 



rj = 2 + r + 



2t 2 



-(2 + 7r + 7r 2 + r 3 ^ 



4(l + r) : 



(11) 



It describes the attraction (repulsion) of a hole at site 
0(1) by the / electron. Furthermore, Hk is a spin-spin 
interaction term in analogy to the cf exchange interac- 
tion, 



H K = 



2jV 2 , 
U-e/ 



tV 2 {2 + r) a ^ 
• S0 + [/e/( l + r) 2 S /- l^&o + s . l 

J i— ±1 



(12) 

with s«/ = (1/2) ^ 0jS ct a cr a ^4/3) where cr a/ 3 are the Pauli 
matrices. The last term H' gives a correction to the 
effective model, 



H' 



2VH 



2,2 



Ue) {l + rf ^ s rMl-™o)-s (l-^)]- (13) 




FIG. 5: Binding energy A B for (a) V = 0.2 and (b) V = 0.1 
with £/ = —3 as a function of the band filling n. The Coulomb 
interaction strengths are U = (crosses), 2 (triangles), 4 (cir- 
cles), and 10 (squares). Filled (empty) symbols correspond 
to the data for n = 1 (n < 1) and empty symbols at n = 1 
represent the values for infinitesimally doped systems. 



The first term of (|13[) implies an antiferromagnetic inter- 
action between the impurity site and site ±1 if there is a 
hole at site 0; on the other hand, the second term gives 
a correction to the Kondo-type interaction, i.e., the first 
term of (fT2|) . and the antiferromagnetic spin exchange 
between the impurity site and site may be reduced. 



1. Binding energy 

Of particular interest is the evolution of the binding 
energy of the impurity-induced bound state upon hole 
doping. We can easily imagine that the binding energy 
is suppressed by hole doping due to the enhancement of 
charge fluctuation. Thus, away from half filling, the ID 
correlated system is metallic and the bound state changes 
from a local singlet to the Kondo singlet. If the bound 
state survives with hole doping, it has a much larger en- 
ergy than the standard Kondo singlet. In Fig. Owe show 
the binding energy Ab as a function of band filling n 
(< 1) at (a) V = 0.2 and (b) V = 0.1 with e f = -3 
for various values of U. Filled (empty) symbols refer 
to the data for n — 1 (n < 1) and empty symbols at 
n = 1 represent the values for infinitesimally doped sys- 
tems (see below). Roughly speaking, Ab is discontin- 
uously reduced at n — 1 and decreases with increasing 
hole doping for all cases except U = 0. We find however 
that Ab remains of the same order of magnitude as in 
the half-filled case even at doping level up to a few per- 
cents. Also, the dependence of Ab on U is weaker for 
higher doping concentrations. 

More precisely, there are two differences in behavior on 
the hybridization strength V. One is that in the vicinity 
of n = 1 the binding energy for V = 0.2 decreases more 
rapidly than that for V — 0.1 despite larger cf exchange 
coupling [Eq.(fT2])]. It must be associated with the at- 
traction between doped holes and the / electron, which is 
described in detail in the next paragraph. The other is 
that the binding energy disappears at lower doping levels 
for small values of V; Ab for V — 0.2 maintains its value 



at n < 0.9 and that for V = 0.1 goes to zero around 
n ~ 0.8 — 0.9. It results from the size of the cf exchange 
coupling J c f, and thus the critical doping concentration 
is highest at U ~ 4 giving a maximal value of J c / . 

For the limit n — > 1 we have extrapolated the finite- 
size binding energy Ab(L) to the thermodynamic limit 
L — ► oo for the four-hole doped system by going up to 
L + 1 = 510. One notices that the value of the binding 
energy in the limit n — > 1 differs from the n = 1 un- 
doped value. It reflects the fact that the binding energy 
of the Kondo singlet in the infinitesimally doped system 
is less than that of the local singlet in the undoped sys- 
tem. The reason being that, when the system is doped 
by a hole, the career tends to move onto site due to 
the attraction from the impurity site [Eq. (|lip ] and thus 
a spin-singlet formation is prevented. The discontinuity 
is higher for V = 0.2 than for V = 0.1 because the at- 
tractive interaction is enhanced by the hybridization V. 
Such a discontinuity of the spin-excitation energy has also 
been found in studies of ladder systems i 17 ' 18 Note that 
in the hole-doped case the V dependence of the binding 
energy is not simple because V enhances two compet- 
ing effects: (i) the attraction between doped holes and 
the / electron and (ii) the cf exchange coupling between 
conduction electrons and the / electron. 



2. Spin-spin correlations 
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FIG. 6: Spin-spin correlation functions (S/ • so) and (S/ ■ si) 
for (a) V — 0.2 and (b) V = 0.1 with £f = —3 as a function 
of the band filling n. The Coulomb interaction strengths are 
U — (crosses), 2 (squares), 4 (triangles), and 10 (circles). 
Filled (empty) symbols correspond to the data for n = 1 
(n < 1) and empty symbols at n = 1 represent the values for 
infinitesimally doped systems. 



Finally, we study the hole-doping dependence of spin- 
spin correlations between the /- and conduction electrons. 
The correlation is expected to be weakened by hole dop- 
ing due to an increase of charge fluctuations. In Fig. O 
we show the spin-spin correlation functions (S/ ■ So) and 
(Sf ■ si) as a function of band filling n (< 1) when (a) 

V = 0.2 and (b) V = 0.1 with Sf = —3 for various 
Coulomb interaction strengths. The properties are fun- 
damentally linked to those of the binding energy as fol- 
lows: (i) correlations are suppressed by hole doping and 
(ii) there exists a discontinuity at n = 1. 

Let us now investigate the DMRG results for the two 

V values. When V — 0.2, all the correlation functions for 
finite U decrease rapidly close to n = 1 and decay slowly 
when n < 0.9. This behavior is quite similar to that 
of the binding energy. It is seen that |(S/ • Si) / (S/ • Sq)| 
decreases with decreasing n. The small value corresponds 
to a rapid decay of (Sf ■ s r ) around r = 0, as seen in 
Fig. 21 e.g., for U — 0.5 and n = 1. It is accompanied by 
a transfer from the local singlet to the Kondo singlet. It 
also suggests a reduction of the RKKY interaction with 
doping. In addition, it is surprising that (S f ■ so) seems 
to be enlarged by hole doping for small values of U (< 2). 
The "exchange hole" around the impurity is as a conse- 
quence of the Pauli principle. When V = 0.1, all the cor- 
relation functions decrease monotonously and go to zero 
around n ~ 0.9, which is accompanied by a vanishing of 
the binding energy. For n < 0.8—0.9, |(S/ • s%) | has small 
negative values for large values of U, which indicates an- 



tiferromagnetic correlations between the / electron and 
the spin at site 1 . It is derived from the first term of (fl3|) 
and was previously suggested in Ref.5. 



V. CONCLUSION 

Using the DMRG method, we have studied a magnetic 
impurity embedded in a correlated electron system which 
is assumed to be the ID Hubbard chain. At half filling, 
we confirm that the binding energy increases exponen- 
tially in the weak-coupling regime. There is a crossover 
from the Kondo singlet to the local singlet. The for- 
mer state involves a wider spread of spin-polarized elec- 
trons around the impurity than the latter one. With in- 
creasing values of U, the binding energy has a maximum 
around U ~ 4 and afterwards decreases inversely propor- 
tional to the Coulomb interaction. Due to the formation 
of a singlet bound state, the spin-spin correlation func- 
tion decays exponentially with distance from the impu- 
rity site for all values of U (> 0). The correlation length 
is quite long when the binding energy is small. It be- 
comes shorter with increasing Coulomb interaction. For 
infinitesimally hole doping, we find a discontinuous re- 
duction of the binding energy and of the spin-spin corre- 
lations from the values at half filling. For further doping, 
the binding energy is reduced but remains of the same 
order of magnitude as in the half-filled case even for dop- 
ing concentration of a few percent. The electron-doped 
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case is not studied here, but we expect qualitatively sini- Acknowledgments 
ilar properties as for hole doping. However, there is no 
discontinuity at half filling. This is so because when an 
electron is added to the half-filled system it is distributed 

uniformly over the ID chain except for site 0. We thank T. Takimoto for useful discussions. 
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